High-throughput identification of functional regulatory SNPs in systemic lupus erythematosus.

Genome-wide association studies implicate multiple loci in risk for systemic lupus erythematosus (SLE), but few contain exonic variants, rendering systematic identification of non-coding variants essential to decoding SLE genetics. We utilized SNP-seq and bioinformatic enrichment to interrogate 2180 single-nucleotide polymorphisms (SNPs) from 87 SLE risk loci for potential binding of transcription factors and related proteins from B cells. 52 SNPs that passed initial screening were tested by electrophoretic mobility shift and luciferase reporter assays. To validate the approach, we studied rs2297550 in detail, finding that the risk allele enhanced binding to the transcription factor Ikaros (IKZF1), thereby modulating expression of IKBKE. Correspondingly, primary cells from genotyped healthy donors bearing the risk allele expressed higher levels of the interferon / NF-κB regulator IKKϵ. Together, these findings define a set of likely functional non-coding lupus risk variants and identify a new regulatory pathway involving rs2297550, Ikaros, and IKKϵ implicated by human genetics in risk for SLE.

Identification of a regulatory pathway governing TRAF1 via an arthritis-associated non-coding variant.

TRAF1/C5 was among the first loci shown to confer risk for inflammatory arthritis in the absence of an associated coding variant, but its genetic mechanism remains undefined. Using Immunochip data from 3,939 patients with juvenile idiopathic arthritis (JIA) and 14,412 control individuals, we identified 132 plausible common non-coding variants, reduced serially by single-nucleotide polymorphism sequencing (SNP-seq), electrophoretic mobility shift, and luciferase studies to the single variant rs7034653 in the third intron of TRAF1. Genetically manipulated experimental cells and primary monocytes from genotyped donors establish that the risk G allele reduces binding of Fos-related antigen 2 (FRA2), encoded by FOSL2, resulting in reduced TRAF1 expression and enhanced tumor necrosis factor (TNF) production. Conditioning on this JIA variant eliminated attributable risk for rheumatoid arthritis, implicating a mechanism shared across the arthritis spectrum. These findings reveal that rs7034653, FRA2, and TRAF1 mediate a pathway through which a non-coding functional variant drives risk of inflammatory arthritis in children and adults.

Molecular pathways identified from single nucleotide polymorphisms demonstrate mechanistic differences in systemic lupus erythematosus patients of Asian and European ancestry.

Systemic lupus erythematosus (SLE) is a multi-organ autoimmune disorder with a prominent genetic component. Individuals of Asian-Ancestry (AsA) disproportionately experience more severe SLE compared to individuals of European-Ancestry (EA), including increased renal involvement and tissue damage. However, the mechanisms underlying elevated severity in the AsA population remain unclear. Here, we utilized available gene expression data and genotype data based on all non-HLA SNP associations in EA and AsA SLE patients detected using the Immunochip genotyping array. We identified 2778 ancestry-specific and 327 trans-ancestry SLE-risk polymorphisms. Genetic associations were examined using connectivity mapping and gene signatures based on predicted biological pathways and were used to interrogate gene expression datasets. SLE-associated pathways in AsA patients included elevated oxidative stress, altered metabolism and mitochondrial dysfunction, whereas SLE-associated pathways in EA patients included a robust interferon response (type I and II) related to enhanced cytosolic nucleic acid sensing and signaling. An independent dataset derived from summary genome-wide association data in an AsA cohort was interrogated and identified similar molecular pathways. Finally, gene expression data from AsA SLE patients corroborated the molecular pathways predicted by SNP associations. Identifying ancestry-related molecular pathways predicted by genetic SLE risk may help to disentangle the population differences in clinical severity that impact AsA and EA individuals with SLE.

Mendelian randomization and pathway analysis demonstrate shared genetic associations between lupus and coronary artery disease.

Coronary artery disease (CAD) is a leading cause of death in patients with systemic lupus erythematosus (SLE). Despite clinical evidence supporting an association between SLE and CAD, pleiotropy-adjusted genetic association studies are limited and focus on only a few common risk loci. Here, we identify a net positive causal estimate of SLE-associated non-HLA SNPs on CAD by traditional Mendelian randomization (MR) approaches. Pathway analysis using SNP-to-gene mapping followed by unsupervised clustering based on protein-protein interactions (PPIs) identifies biological networks composed of positive and negative causal sets of genes. In addition, we confirm the casual effects of specific SNP-to-gene modules on CAD using only SNP mapping to each PPI-defined functional gene set as instrumental variables. This PPI-based MR approach elucidates various molecular pathways with causal implications between SLE and CAD and identifies biological pathways likely causative of both pathologies, revealing known and novel therapeutic interventions for managing CAD in SLE.

Gut dysbiosis and the clinical spectrum in anti-Ro positive mothers of children with neonatal lupus.

Anti-SSA/Ro antibodies, while strongly linked to fetal cardiac injury and neonatal rash, can associate with a spectrum of disease in the mother, ranging from completely asymptomatic to overt Systemic Lupus Erythematosus (SLE) or Sjögren's Syndrome (SS). This study was initiated to test the hypothesis that the microbiome, influenced in part by genetics, contributes to disease state. The stool microbiome of healthy controls (HC) was compared to that of anti-SSA/Ro positive women whose children had neonatal lupus. At the time of sampling, these women were either asymptomatic (Asym), had minor rheumatic symptoms or signs considered as an undifferentiated autoimmune syndrome (UAS), or were diagnosed with SLE or SS. Differences in microbial relative abundances among these three groups were tested assuming an ordering in clinical severity (HC

Nucleic Acid-Sensing and Interferon-Inducible Pathways Show Differential Methylation in MZ Twins Discordant for Lupus and Overexpression in Independent Lupus Samples: Implications for Pathogenic Mechanism and Drug Targeting.

Systemic lupus erythematosus (SLE) is a chronic, multisystem, autoimmune inflammatory disease with genomic and non-genomic contributions to risk. We hypothesize that epigenetic factors are a significant contributor to SLE risk and may be informative for identifying pathogenic mechanisms and therapeutic targets. To test this hypothesis while controlling for genetic background, we performed an epigenome-wide analysis of DNA methylation in genomic DNA from whole blood in three pairs of female monozygotic (MZ) twins of European ancestry, discordant for SLE. Results were replicated on the same array in four cell types from a set of four Danish female MZ twin pairs discordant for SLE. Genes implicated by the epigenetic analyses were then evaluated in 10 independent SLE gene expression datasets from the Gene Expression Omnibus (GEO). There were 59 differentially methylated loci between unaffected and affected MZ twins in whole blood, including 11 novel loci. All but two of these loci were hypomethylated in the SLE twins relative to the unaffected twins. The genes harboring these hypomethylated loci exhibited increased expression in multiple independent datasets of SLE patients. This pattern was largely consistent regardless of disease activity, cell type, or renal tissue type. The genes proximal to CpGs exhibiting differential methylation (DM) in the SLE-discordant MZ twins and exhibiting differential expression (DE) in independent SLE GEO cohorts (DM-DE genes) clustered into two pathways: the nucleic acid-sensing pathway and the type I interferon pathway. The DM-DE genes were also informatically queried for potential gene-drug interactions, yielding a list of 41 drugs including a known SLE therapy. The DM-DE genes delineate two important biologic pathways that are not only reflective of the heterogeneity of SLE but may also correlate with distinct IFN responses that depend on the source, type, and location of nucleic acid molecules and the activated receptors in individual patients. Cell- and tissue-specific analyses will be critical to the understanding of genetic factors dysregulating the nucleic acid-sensing and IFN pathways and whether these factors could be appropriate targets for therapeutic intervention.

Combined genetic analysis of juvenile idiopathic arthritis clinical subtypes identifies novel risk loci, target genes and key regulatory mechanisms.

OBJECTIVES: Juvenile idiopathic arthritis (JIA) is the most prevalent form of juvenile rheumatic disease. Our understanding of the genetic risk factors for this disease is limited due to low disease prevalence and extensive clinical heterogeneity. The objective of this research is to identify novel JIA susceptibility variants and link these variants to target genes, which is essential to facilitate the translation of genetic discoveries to clinical benefit. METHODS: We performed a genome-wide association study (GWAS) in 3305 patients and 9196 healthy controls, and used a Bayesian model selection approach to systematically investigate specificity and sharing of associated loci across JIA clinical subtypes. Suggestive signals were followed-up for meta-analysis with a previous GWAS (2751 cases/15 886 controls). We tested for enrichment of association signals in a broad range of functional annotations, and integrated statistical fine-mapping and experimental data to identify target genes. RESULTS: Our analysis provides evidence to support joint analysis of all JIA subtypes with the identification of five novel significant loci. Fine-mapping nominated causal single nucleotide polymorphisms with posterior inclusion probabilities ≥50% in five JIA loci. Enrichment analysis identified RELA and EBF1 as key transcription factors contributing to disease risk. Our integrative approach provided compelling evidence to prioritise target genes at six loci, highlighting mechanistic insights for the disease biology and IL6ST as a potential drug target. CONCLUSIONS: In a large JIA GWAS, we identify five novel risk loci and describe potential function of JIA association signals that will be informative for future experimental works and therapeutic strategies.

Platelet Glycoprotein Ib α-Chain as a Putative Therapeutic Target for Juvenile Idiopathic Arthritis: A Mendelian Randomization Study.

OBJECTIVE: To ascertain the role of platelet glycoprotein Ib α-chain (GPIbα) plasma protein levels in cardiovascular, autoimmune, and autoinflammatory diseases and whether its effects are mediated by platelet count. METHODS: We performed a two-sample Mendelian randomization (MR) study, using both a cis-acting protein quantitative trait locus (cis-pQTL) and trans-pQTL near the GP1BA and BRAP genes as instruments. To assess if platelet count mediated the effect, we then performed a two-step MR study. Putative associations (GPIbα/platelet count/disease) detected by MR analyses were subsequently assessed using multiple-trait colocalization analyses. RESULTS: After correction for multiple testing (Bonferroni-corrected threshold P ≤ 2 × 10-3 ), GPIbα, instrumented by either cis-pQTL or trans-pQTL, was causally implicated with an increased risk of oligoarticular and rheumatoid factor (RF)-negative polyarticular juvenile idiopathic arthritis (JIA). These effects of GPIbα appeared to be mediated by platelet count and were supported by strong evidence of colocalization (probability of all 3 traits sharing a common causal variant ≥0.80). GPIbα instrumented by cis-pQTL did not appear to affect cardiovascular risk, although the GPIbα trans-pQTL was associated with an increased risk of cardiovascular diseases and autoimmune diseases but a decreased risk of autoinflammatory diseases, suggesting that this trans-acting instrument operates through other pathways. CONCLUSION: The role of platelets in thrombosis is well-established; however, our findings provide some novel genetic evidence that platelets may be causally implicated in the development of oligoarticular and RF-negative polyarticular JIA, and indicate that GPIbα may serve as a putative therapeutic target for these JIA subtypes.

Analysis of Trans-Ancestral SLE Risk Loci Identifies Unique Biologic Networks and Drug Targets in African and European Ancestries.

Systemic lupus erythematosus (SLE) is a multi-organ autoimmune disorder with a prominent genetic component. Individuals of African ancestry (AA) experience the disease more severely and with an increased co-morbidity burden compared to European ancestry (EA) populations. We hypothesize that the disparities in disease prevalence, activity, and response to standard medications between AA and EA populations is partially conferred by genomic influences on biological pathways. To address this, we applied a comprehensive approach to identify all genes predicted from SNP-associated risk loci detected with the Immunochip. By combining genes predicted via eQTL analysis, as well as those predicted from base-pair changes in intergenic enhancer sites, coding-region variants, and SNP-gene proximity, we were able to identify 1,731 potential ancestry-specific and trans-ancestry genetic drivers of SLE. Gene associations were linked to upstream and downstream regulators using connectivity mapping, and predicted biological pathways were mined for candidate drug targets. Examination of trans-ancestral pathways reflect the well-defined role for interferons in SLE and revealed pathways associated with tissue repair and remodeling. EA-dominant genetic drivers were more often associated with innate immune and myeloid cell function pathways, whereas AA-dominant pathways mirror clinical findings in AA subjects, suggesting disease progression is driven by aberrant B cell activity accompanied by ER stress and metabolic dysfunction. Finally, potential ancestry-specific and non-specific drug candidates were identified. The integration of all SLE SNP-predicted genes into functional pathways revealed critical molecular pathways representative of each population, underscoring the influence of ancestry on disease mechanism and also providing key insight for therapeutic selection.

Effects of Intensive Systolic Blood Pressure Control on All-Cause Hospitalizations.

Intensive blood pressure control decreases the rate of cardiovascular events by >25% compared with standard blood pressure control. We sought to determine whether the decrease in cardiovascular events seen with intensive blood pressure control is associated with an increased rate of other causes of hospitalization. This is a post hoc analysis of SPRINT (Systolic Blood Pressure Intervention Trial) in 9361 adult participants with hypertension and elevated cardiovascular risk. Participants were randomly assigned to an intensive or standard systolic blood pressure goal (<120 or <140 mm Hg, respectively). The primary outcome was hospitalization rates per 100 person-years for hospitalizations not associated with SPRINT primary events. After excluding hospitalizations linked to SPRINT primary events, there were 4678 participants with a rate of 19.70 hospitalizations per 100 person-years, compared with 4683 participants with a rate of 19.65 (P=0.37). Equivalence testing shows that these hospitalization rates were statistically equivalent at the P=0.05 level. Of those with hospitalizations, >1 hospitalization was seen in 38.8% of intensive arm participants and 41.9% of standard arm participants (P=0.08). The mean cumulative count of nonprimary event hospitalizations was comparable between the two arms. The most common causes of hospitalization were cardiovascular (23.6%) followed by injuries, including bone and joint therapeutic procedures (15.7%), infections (12.0%), and nervous systems disorders (10.7%). No categories of hospitalization were statistically more common in the intensive arm compared with the standard arm. Thus, the decrease in cardiovascular events seen with intensive blood pressure control is not associated with an increased rate of other causes of hospitalization. Registration- URL: https://www.clinicaltrials.gov; Unique identifier: NCT01206062.

Association of SLCO1B1 *14 Allele with Poor Response to Methotrexate in Juvenile Idiopathic Arthritis Patients.

OBJECTIVE: Variants in the SLCO1B1 gene, encoding a hepatic methotrexate (MTX) transporter, affect clearance of high-dose MTX. We tested whether in the *14 and *15 alleles of SLCO1B1 influenced the response to low-dose MTX in juvenile idiopathic arthritis (JIA) patients. METHODS: The study included 310 JIA patients genotyped for three single nucleotide polymorphisms (SNPs) in SLCO1B1 (rs4149056, rs2306283, and rs11045819). A patient's SLCO1B1 diplotype was determined by combining the SNPs into the *1a, *1b, *4, *5, *14, and *15 alleles. Number of active joints at follow-up (visit closest to 6 months of treatment and prior to starting a tumor necrosis factor inhibitor) was used as the dependent variable in a negative binomial regression model that included active joint count at baseline as a covariate. RESULTS: The SLCO1B1*14 allele was associated with less response to MTX (P = 0.024) and the *15 allele was not associated with response to MTX (P = 0.392). CONCLUSION: SLCO1B1 alleles may be associated with poor response to MTX in JIA patients. The *14 allele has been associated with fast clearance (low exposure) after high-dose MTX in patients with leukemia. Thus, the SLCO1B1 gene may be informative for precision dosing of MTX in JIA patients. Patients carrying the *14 allele may require a higher dose than noncarriers to achieve a similar response to MTX.

Trans-Ethnic Mapping of BANK1 Identifies Two Independent SLE-Risk Linkage Groups Enriched for Co-Transcriptional Splicing Marks.

BANK1 is a susceptibility gene for several systemic autoimmune diseases in several populations. Using the genome-wide association study (GWAS) data from Europeans (EUR) and African Americans (AA), we performed an extensive fine mapping of ankyrin repeats 1 (BANK1). To increase the SNP density, we used imputation followed by univariate and conditional analysis, combined with a haplotypic and expression quantitative trait locus (eQTL) analysis. The data from Europeans showed that the associated region was restricted to a minimal and dependent set of SNPs covering introns two and three, and exon two. In AA, the signal found in the Europeans was split into two independent effects. All of the major risk associated SNPs were eQTLs, and the risks were associated with an increased BANK1 gene expression. Functional annotation analysis revealed the enrichment of repressive B cell epigenomic marks (EZH2 and H3K27me3) and a strong enrichment of splice junctions. Furthermore, one eQTL located in intron two, rs13106926, was found within the binding site for RUNX3, a transcriptional activator. These results connect the local genome topography, chromatin structure, and the regulatory landscape of BANK1 with co-transcriptional splicing of exon two. Our data defines a minimal set of risk associated eQTLs predicted to be involved in the expression of BANK1 modulated through epigenetic regulation and splicing. These findings allow us to suggest that the increased expression of BANK1 will have an impact on B-cell mediated disease pathways.

A plausibly causal functional lupus-associated risk variant in the STAT1-STAT4 locus.

Systemic lupus erythematosus (SLE or lupus) (OMIM: 152700) is a chronic autoimmune disease with debilitating inflammation that affects multiple organ systems. The STAT1-STAT4 locus is one of the first and most highly replicated genetic loci associated with lupus risk. We performed a fine-mapping study to identify plausible causal variants within the STAT1-STAT4 locus associated with increased lupus disease risk. Using complementary frequentist and Bayesian approaches in trans-ancestral Discovery and Replication cohorts, we found one variant whose association with lupus risk is supported across ancestries in both the Discovery and Replication cohorts: rs11889341. In B cell lines from patients with lupus and healthy controls, the lupus risk allele of rs11889341 was associated with increased STAT1 expression. We demonstrated that the transcription factor HMGA1, a member of the HMG transcription factor family with an AT-hook DNA-binding domain, has enriched binding to the risk allele compared with the non-risk allele of rs11889341. We identified a genotype-dependent repressive element in the DNA within the intron of STAT4 surrounding rs11889341. Consistent with expression quantitative trait locus (eQTL) analysis, the lupus risk allele of rs11889341 decreased the activity of this putative repressor. Altogether, we present a plausible molecular mechanism for increased lupus risk at the STAT1-STAT4 locus in which the risk allele of rs11889341, the most probable causal variant, leads to elevated STAT1 expression in B cells due to decreased repressor activity mediated by increased binding of HMGA1.

Brief Report: The Genetic Profile of Rheumatoid Factor-Positive Polyarticular Juvenile Idiopathic Arthritis Resembles That of Adult Rheumatoid Arthritis.

OBJECTIVE: Juvenile idiopathic arthritis (JIA) comprises 7 heterogeneous categories of chronic childhood arthritides. Approximately 5% of children with JIA have rheumatoid factor (RF)-positive arthritis, which phenotypically resembles adult rheumatoid arthritis (RA). Our objective was to compare and contrast the genetics of RF-positive polyarticular JIA with those of RA and selected other JIA categories, to more fully understand the pathophysiologic relationships of inflammatory arthropathies. METHODS: Patients with RF-positive polyarticular JIA (n = 340) and controls (n = 14,412) were genotyped using the Immunochip array. Single-nucleotide polymorphisms were tested for association using a logistic regression model adjusting for admixture proportions. We calculated weighted genetic risk scores (wGRS) of reported RA and JIA risk loci, and we compared the ability of these wGRS to predict RF-positive polyarticular JIA. RESULTS: As expected, the HLA region was strongly associated with RF-positive polyarticular JIA (P = 5.51 × 10-31 ). Nineteen of 44 RA risk loci and 6 of 27 oligoarticular/RF-negative polyarticular JIA risk loci were associated with RF-positive polyarticular JIA (P < 0.05). The RA wGRS predicted RF-positive polyarticular JIA (area under the curve [AUC] 0.71) better than did the oligoarticular/RF-negative polyarticular JIA wGRS (AUC 0.59). The genetic profile of patients with RF-positive polyarticular JIA was more similar to that of RA patients with age at onset 16-29 years than to that of RA patients with age at onset ≥70 years. CONCLUSION: RF-positive polyarticular JIA is genetically more similar to adult RA than to the most common JIA categories and thus appears to be a childhood-onset presentation of autoantibody-positive RA. These findings suggest common disease mechanisms, which could lead to novel therapeutic targets and shared treatment strategies.

Association of Natural Killer Cell Ligand Polymorphism HLA-C Asn80Lys With the Development of Anti-SSA/Ro-Associated Congenital Heart Block.

OBJECTIVE: Fetal exposure to maternal anti-SSA/Ro antibodies is necessary but not sufficient for the development of autoimmune congenital heart block (CHB), suggesting that other factors, such as fetal genetic predisposition, are important. Given the previously described association between major histocompatibility complex alleles and CHB risk, we undertook the present study to test the hypothesis that a variant form of HLA-C Asn80Lys, which binds with high affinity to an inhibitory killer cell immunoglobulin-like receptor (KIR) and thus renders natural killer (NK) cells incapable of restricting inflammation, contributes to the development of CHB. METHODS: Members of 192 pedigrees in the US and Europe (194 cases of CHB, 91 unaffected siblings, 152 fathers, 167 mothers) and 1,073 out-of-study controls were genotyped on the Immunochip single-nucleotide polymorphism microarray. Imputation was used to identify associations at HLA-C Asn80Lys (Asn, C1; Lys, C2) and KIR. Tests for association were performed using logistic regression. McNemar's test and the pedigree disequilibrium test (PDT) were used for matched analyses between affected and unaffected children. RESULTS: Compared with out-of-study controls of the same sex, the C2 allele was less frequent in the mothers (odds ratio [OR] 0.63, P = 0.0014) and more frequent in the fathers (OR 1.40, P = 0.0123), yielding a significant sex-by-C2 interaction (P = 0.0002). The C2 allele was more frequent in affected siblings than in unaffected siblings (OR 3.67, P = 0.0025), which was consistent with the PDT results (P = 0.016); these results were observed in both sexes and across the US and European cohorts. There was no difference in the frequency of the inhibitory KIR genotype (KIR AA) between affected and unaffected children (P = 0.55). CONCLUSION: These data establish C2 as a novel genetic risk factor associated with CHB. This observation supports a model in which fetuses with C2 ligand expression and maternal anti-SSA/Ro positivity may have impaired NK cell surveillance, resulting in unchecked cardiac inflammation and scarring.

Transancestral mapping and genetic load in systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) is an autoimmune disease with marked gender and ethnic disparities. We report a large transancestral association study of SLE using Immunochip genotype data from 27,574 individuals of European (EA), African (AA) and Hispanic Amerindian (HA) ancestry. We identify 58 distinct non-HLA regions in EA, 9 in AA and 16 in HA (∼50% of these regions have multiple independent associations); these include 24 novel SLE regions (P<5 × 10-8), refined association signals in established regions, extended associations to additional ancestries, and a disentangled complex HLA multigenic effect. The risk allele count (genetic load) exhibits an accelerating pattern of SLE risk, leading us to posit a cumulative hit hypothesis for autoimmune disease. Comparing results across the three ancestries identifies both ancestry-dependent and ancestry-independent contributions to SLE risk. Our results are consistent with the unique and complex histories of the populations sampled, and collectively help clarify the genetic architecture and ethnic disparities in SLE.

Genome-Wide Association Meta-Analysis Reveals Novel Juvenile Idiopathic Arthritis Susceptibility Loci.

OBJECTIVE: Juvenile idiopathic arthritis (JIA) is the most common childhood rheumatic disease and has a strong genomic component. To date, JIA genetic association studies have had limited sample sizes, used heterogeneous patient populations, or included only candidate regions. The aim of this study was to identify new associations between JIA patients with oligoarticular disease and those with IgM rheumatoid factor (RF)-negative polyarticular disease, which are clinically similar and the most prevalent JIA disease subtypes. METHODS: Three cohorts comprising 2,751 patients with oligoarticular or RF-negative polyarticular JIA were genotyped using the Affymetrix Genome-Wide SNP Array 6.0 or the Illumina HumanCoreExome-12+ Array. Overall, 15,886 local and out-of-study controls, typed on these platforms or the Illumina HumanOmni2.5, were used for association analyses. High-quality single-nucleotide polymorphisms (SNPs) were used for imputation to 1000 Genomes prior to SNP association analysis. RESULTS: Meta-analysis showed evidence of association (P < 1 × 10-6 ) at 9 regions: PRR9_LOR (P = 5.12 × 10-8 ), ILDR1_CD86 (P = 6.73 × 10-8 ), WDFY4 (P = 1.79 × 10-7 ), PTH1R (P = 1.87 × 10-7 ), RNF215 (P = 3.09 × 10-7 ), AHI1_LINC00271 (P = 3.48 × 10-7 ), JAK1 (P = 4.18 × 10-7 ), LINC00951 (P = 5.80 × 10-7 ), and HBP1 (P = 7.29 × 10-7 ). Of these, PRR9_LOR, ILDR1_CD86, RNF215, LINC00951, and HBP1 were shown, for the first time, to be autoimmune disease susceptibility loci. Furthermore, associated SNPs included cis expression quantitative trait loci for WDFY4, CCDC12, MTP18, SF3A1, AHI1, COG5, HBP1, and GPR22. CONCLUSION: This study provides evidence of both unique JIA risk loci and risk loci overlapping between JIA and other autoimmune diseases. These newly associated SNPs are shown to influence gene expression, and their bounding regions tie into molecular pathways of immunologic relevance. Thus, they likely represent regions that contribute to the pathology of oligoarticular JIA and RF-negative polyarticular JIA.

Lupus Risk Variant Increases pSTAT1 Binding and Decreases ETS1 Expression.

Genetic variants at chromosomal region 11q23.3, near the gene ETS1, have been associated with systemic lupus erythematosus (SLE), or lupus, in independent cohorts of Asian ancestry. Several recent studies have implicated ETS1 as a critical driver of immune cell function and differentiation, and mice deficient in ETS1 develop an SLE-like autoimmunity. We performed a fine-mapping study of 14,551 subjects from multi-ancestral cohorts by starting with genotyped variants and imputing to all common variants spanning ETS1. By constructing genetic models via frequentist and Bayesian association methods, we identified 16 variants that are statistically likely to be causal. We functionally assessed each of these variants on the basis of their likelihood of affecting transcription factor binding, miRNA binding, or chromatin state. Of the four variants that we experimentally examined, only rs6590330 differentially binds lysate from B cells. Using mass spectrometry, we found more binding of the transcription factor signal transducer and activator of transcription 1 (STAT1) to DNA near the risk allele of rs6590330 than near the non-risk allele. Immunoblot analysis and chromatin immunoprecipitation of pSTAT1 in B cells heterozygous for rs6590330 confirmed that the risk allele increased binding to the active form of STAT1. Analysis with expression quantitative trait loci indicated that the risk allele of rs6590330 is associated with decreased ETS1 expression in Han Chinese, but not other ancestral cohorts. We propose a model in which the risk allele of rs6590330 is associated with decreased ETS1 expression and increases SLE risk by enhancing the binding of pSTAT1.